RESPONSE OF A MODEL OF CO OXIDATION WITH CO DESORPTION AND DIFFUSION TO 

A PERIODIC EXTERNAL CO PRESSURE 



o 



O 
(N 



C/3 



G. M.Buendia and E. Machado 
Physics Department, Universidad Simon Bolivar, 
Apartado 89000, Caracas 1080, Venezuela 



O ■ P. A. Rikvold 



Center for Materials Research and Technology, 
School of Computational Science, and Department of Physics, 
^ , Florida State University, Tallahassee, Florida 32306-4-120, USA 



Abstract 

We present a study of the dynamical behavior of a Ziff-Gulari-Barshad model with CO desorption and lateral 
diffusion. Depending on the values of the desorption and diffusion parameters, the system presents a discontinuous 
phase transition between low and high CO coverage phases. We calculate several points on the coexistence curve 
between these phases. Inclusion of the diffusion term produces a significant increase in the CO2 production rate. We 
further applied a square- wave periodic pressure variation of the partial CO pressure with parameters that can be tuned 
to modify the catalytic activity. Contrary to the diffusion-free case, this driven system does not present a further 
enhancement of the catalytic activity, beyond the increase induced by the diffusion under constant CO pressure. 

I. INTRODUCTION 



. Phase transitions and critical phenomena in open noncquilibrium statistical systems have attracted a great deal of 
O ' interest in recent decades, due to their wide applications in many branches of physics, chemistry, biology, economics, 
, and even sociology T|. Among these, surface reaction models have become the archetype to study out-of-equilibrium 

critical phenomena, and they have been intensely analyzed with the purpose of designing more efficient catalytic 
processes 0. The Ziff, Gulari, Barshad (ZGB) model [g with desorption (ZGB-k) HlllEla 

describes some kinetic 

K*" ' aspects of the reaction CO-f O2 — > CO2 on a catalytic surface in terms of two parameters: the partial pressure of CO, 
' y, which represents the probability per unit of time that the next molecule arriving at the surface is a CO, and the 
desorption parameter, fc, which gives the probability per unit of time that a site occupied by a CO molecule is vacated. 
Below a critical CO desorption rate, the model exhibits a first-order phase transition between low and high coverage 
. phases 0, 0, S 13 • A study of the metastable lifetimes associated with the transition between these phases indicates 
\^ ' that the decay mechanism is very similar to the decay of metastable phases associated with equilibrium first-order 
phase transitions and can be described by the classic Kolmorogov-Johnson-Mehl-Avrami (KJMA) theory 9] of phase 
transformations by nucleation and growth 0. Near the coexistence line between the active and the CO poisoned 
phase, the decay times of the metastable state are different if the ZGB-k model is driven into the active phase from 
H the CO poisoned state or vice versa 0- This asymmetry of the lifetimes, together with studies that indicate that 
' it is possible to increase the efficiency of catalytic reactions by forcing the system by a periodic external pressure 
^5 that switches back and forth around the discontinuous transition inspired us to subject the system to a 

Q square-wave periodic CO pressure that stays for a time td in the high-production region and tp in the low-production 
^ , region. We showed that tp and td can be tuned to significantly enhance the time-averaged catalytic activity of the 
•• • system @. 

. ^ [ Many experiments indicate that at intermediate and high temperatures adsorbed CO is fairly mobile, while oxygen 
is still relatively immobile '^,'10]. In this work we aim at a more realistic process by incorporating CO diffusion into 
?H the ZGB-k model. The model that includes CO desorption and diffusion will be denoted the ZGB-(k,d) model. We 
investigate how the diffusion process affects the first-order kinetic phase transition and the response of the system to 
a periodic variation of the CO pressure. 

The rest of this paper is organized as follows: in Sec. II we define the model and describe the Monte Carlo 
simulation techniques used. In Sec. Ill we present and discuss the results obtained when subjecting the model to a 
periodic external CO pressure. Finally, we present our conclusions in Sec. IV. 
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II. MODEL AND SIMULATIONS 



The ZGB-(k,d) model is simulated on a square lattice of linear size L that represents the catalyst surface. A Monte 
Carlo simulation generates a sequence of trials: CO diffusion with probability per unit of time d, and adsorption 
or desorption with (1 — d). In nature, the diffusion process occurs more frequently than the adsorption/desorption 
process, so that d e (1/2,1) The chemical processes include CO or O2 adsorption with probability per unit 

of time 1 — /c, and CO desorption with probability per unit of time k. For the adsorption, a CO or O2 molecule is 
selected with probability per unit of time y and 1 — y respectively. The algorithm works in the following way. A site i 
is selected at random. In the case of CO diffusion, if i is empty or occupied by an O atom, the trial ends. Otherwise, 
a nearest neighbor of i is selected at random; if it is vacant, the CO molecule exchanges with it, and reacts if any 
of its nearest-neighbors is an O atom, producing a CO2 molecule and leaving two empty sites on the surface. In the 
case of desorption, if i is occupied by a CO, the site i is vacated and the trial ends; if not, the trial also ends. In the 
case of adsorption, if a CO molecule is selected it can be adsorbed at the empty site i if none of its nearest neighbors 
are occupied by an O atom. Otherwise, one of the occupied O neighbors is selected at random and removed from the 
surface, leaving i and the selected neighbor vacant. O2 molecules require a nearest-neighbor pair of vacant sites to 
adsorb. Once the O2 molecule is adsorbed, it is dissociated into two O atoms. If an O atom is located next to a site 
filled with a CO molecule, they react to form a CO2 molecule that escapes, leaving two sites vacant. This process 
simulates the CO -I- O ^ CO2 surface reaction. Since the diffusion is much faster than the other chemical processes, 
we will assume that the diffusion does not slow down the time scale for adsorption/desorption, which is the physically 
relevant time scale for the phase transformation Thus, we define the time unit as one Monte Carlo Step per Site, 
MCSS, in which the system performs, on average, one adsorption/desorption attempt per site. 



III. RESULTS 



Our simulations are reahzed on a square lattice of 100x100 sites, assuming periodic boundary conditions. We 
calculate the coverages 6*00 and 6*0, defined as the fraction of surface sites occupied by CO and O, respectively, and 
i?C02 defined as the rate of production of CO2. In Fig. ^ we show, for a particular value of k and d, the dependence 
of the average values of the coverages and the production rate: (^co), (^o), and {Rco2)i respectively, for an external 
constant CO pressure y. There are two inactive regions, y < yi (yi appears to be fairly independent of k and <i |3])j 
and y > y2{k,d). These correspond to the cases in which the surface is saturated with O and CO, respectively. In 
between, for yi < y < 2/2(^7 d), there is an active window where the system produces CO2. Notice that the maximum 
value of {Rco^) is reached as y2{k,d) is approached from the low-6'co phase. Fig. ^b) shows that the effect of the 
CO diffusion is to widen the range of values of y for which there is catalytic activity, producing a significant increase 
of the maximum CO2 production rate. In Fig.|21we plot the coexistence point, 2/2, as a function of k for several values 
of d. The effect of the diffusion is to shift the transition to higher values of ?/2, while maintaining approximately the 
same dependence of j/2 on k that is observed at d = 0. 

Next, we perturb the system by applying an oscillating pressure y{t) that in a period T = td + tp takes the values , 

{ yi during the time interval tj, 
^ \ yh during the time interval tp , 

located at both sides of the transition point, i.e., yi < y2{k) < yh (see Fig.|2Ia)). 

We found that, for each choice of yi and yh, we can vary the productivity of the system by tuning td and tp, the 
times that the driving force spends in the high and low coverage regions, respectively. In Fig.|3|[b) it can be seen that 
the response to the periodic pressure shown in Fig. |3fa), the CO2 production rate i?c02J also exhibits an oscillatory 
behavior. We define the period-averaged value of i?c02 > 

r ^ ^ J Rco,{t)dt, (2) 

as the dynamic order parameter. For the parameters used in Fig. |31 (selected after an exhaustive search looking 
for the maximum productivity) the long-time average of r, (r) is almost the same as the maximum average CO2 
production rate for constant y. In contrast, the ZGB-k model shows a marked enhancement of the catalytic activity 
when subjected to this type of periodic perturbation |^. We think that this significant difference may possibly be due 
to the fact that the diffusion term modifies the decay mechanisms of the system near the transition. Previous studies 
in a lattice-gas model near an equilibrium first-order phase transition show that at high diffusion rates the nucleation 
no longer can be described by the classical KJMA theory The fact that applying a periodic CO pressure to the 
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system with diffusion, does not seem to increase tfie CO2 production further, may indicate that the optimal mixing 
of CO and O on the surface was already reached by the inclusion of diffusion. 

Depending on the values of td and tp, the system behavior defines two regions: a productive one with r ^ 0, and 
a non-productive one with r w 0, separated by a coexistence curve. In Fig. 2|we show a density plot of the average 
productivity (r) in terms of and tp for a particular choice of yi and yh- Notice that the highest productivity is 
reached close to the coexistence curve. The plot suggests, as is also the case in the ZGB-k model 'Gj, the existence 
of a dynamic phase transition (DPT), where the period-averaged order parameter r takes the value r « or r 7^ 
depending on whether the system is in a nonproductive state or in a productive one. A more detailed study of the 
nature of the transition, based on a finite-size scaling analysis of the fluctuations of the order parameter, is in progress 
and will be published elsewhere. 

IV. CONCLUSIONS 

In this paper we have studied the ZGB model with CO desorption and diffusion, the ZGB-(k,d) model, and its 
dynamic response to periodic variations of the partial CO pressure. We found that, like the model without diffusion 
(ZGB-k), this system exhibits a non-equilibrium discontinuous phase transition between low and high CO coverage 
phases. The transition occurs at a /c-dependent value of the CO pressure 2/2- The effect of the diffusion term is to shift 
the transition to a higher pressure, such that ?/2(fc, d) > y2{k, 0), and to produce a significant increase of the maximum 
CO2 production rate. We further subjected the system to a square-wave periodic CO pressure that switches back 
and forth around the coexistence value that separates the low and high CO-coverage phases, staying for a time td 
in the high-production region and for a time tp in the low production one. We found that tp and td can be tuned 
to vary the catalytic activity of the system. However, contrary to the case of the ZGB-k model, the increase in the 
production of CO2 obtained in this way is negligible. We believe that this result indicates that the decay mechanisms 
of the metastable phase may be strongly affected by the diffusion process. We also found indications that this driven 
nonequilibrium system could have a dynamic phase transition between a dynamic phase of high CO2 production and a 
nonproductive phase. A more detailed study that includes the application of finite-size scaling techniques is necessary 
to establish the full details of the transition. 
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FIG. 1: (a) Average values of the CO and O coverages, ^co and 9o, shown as functions of the stationary apphed CO pressure 
y, for L — 100 with k = 0.01 and d = 0.9. A continuous, nonequilibriurn phase transition occurs at j/i, and a discontinuous one 
at 2/2 (A:, d). (b) Comparison between the average CO2 production rate (iicog) at d = 0.9 and d = 0, calculated for the same 
values of L and as in (a). Both ?/2 and the maximum value of (Jico2) increase significantly with d. 
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FIG. 2: Some points on the coexistence curve for several values of the diffusion probability per unit of time d. This is analogous 
to the pressure vs temperature phase diagram for a fluid in equilibrium. The continuous dashed lines represent a quadratic fit 
to the data. 
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FIG. 3: (a) Applied periodic pressure of CO, y{t), which takes the values yi = 0.54 and yn ~ 0.59 during the time intervals 
td = 10 and tp = 30, respectively, (b) Response of the CO2 production rate to the applied pressure given in (a) for k = 0.01 
and d = 0.9. The dotted line marked (r) indicates the long-time average of the period-averaged CO2 production rate r. In 
contrast to the case for d — 0, this value cannot be distinguished from the meiximum average CO2 production rate for constant 
y ~ 2/2 (fc, d), shown in Fig.0 Time is measured in units of Monte Carlo steps per site. 
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FIG. 4: Long-time average (r) of the period-averaged CO2 production rate, shown as a density plot vs td and tp for yi = 0.54, 
yn = 0.59, k = 0.01, d = 0.9 and L = 100. 
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